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Rational  design  methods  based  on  a  statistical 
description  of  both  the  ship  load  and  strength  have  been 
proposed  by  several  authors  in  lieu  of  the  traditional 
bending  moment  calculations.   In  general,  the  effect  of 
horizontal  moments  on  the  primary  strength  has  been 
neglected.   It  is  shown  in  this  paper,  as  several  others 
have  found,  that  the  horizontal  moments  can  be  very 
significant.   Two  methods  are  described  to  obtain  the 
magnitude  of  the  equivalent  combined  moment. 

An  example  application  is  shown  for  a  large  tanker 
by  use  of  the  MIT  seakeeping  program.   Wave  moment  trends 
with  speed,  sea  state  and  wave  heading  are  noted.   Measured 
and  predicted  values  of  the  vertical  bending  stress  are 
compared.   Finally,  the  parameters  of  the  weibull  distrib- 
ution to  describe  long-term  wave  moment  amplitudes  are 
derived. 
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NOTATION 

L  =  ship's  length 

H  =  wave  height 

h  =  significant  wave  height 

M^  =  vertical  bending  moment 

M     =  horizontal  bending  moment 
H 

M     =  equivalent  combined  bending  moment 
k..    =  vertical  section  modulus 
k„    =  horizontal  section  modulus 

n 
K     =  section  modulus  ratio 

E  =  mean  square  value  of  the  amplitude  of  the  response 

s,k  -  parameters  of  the  Weibull  distribution 

m  =  mean  value  of  the  response  amplitude 

V  =  amplitude  of  the  vertical  moment  component 

H     =  amplitude  of  the  horizontal  moment  component 

V  =  ship's  speed 

L  =  effective  wave  length 

RAO  =  response  amplitude  operator 

RMS  =  root  mean  square 

u  =  circular  frequency 

0      =  phase  angle  of  the  horizontal  moment  with  respect  to  a 
wave  with  its  crest  amidship 

e     =  phase  angle  of  the  vertical  moment  with  respect  to 
a  wave  with  its  crest  amidship 

<j>     =  phase  angle  of  the  equivalent  combined  moment 

A.    =  correlation  coefficient  of  the  vertical  and  horizontal 
bending  moment 

vi 


aT  =  total  or  combined  stress 

av  =  vertical  bending  stress 

a„  =  horizontal  bending  stress 

H 

X  =  wave  length 

S^tw)  =  spectral  density  of  the  response 

H(cj)  =  RAO  of  the  ship 

S(w)  =  sea  spectral  density 
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Chapter  I 

INTRODUCTION 

The  traditional  method  for  determining  the  load  for 
primary  strength  calculations  consists  of  poising  the  ship 
on  a  trochoidal  wave  of  length  equal  to  that  of  the  ship  and 
wave  height  equal  to  one-twentieth  (1/20)  of  the  length. 
The  static  bending  moment  is  then  determined  by  the  load- 
shear-moment  relations  from  beam  theory.   There  have  been 

n  fi 
various  proposals  to  modify  the  wave  height  (e.g.,  H  =  0.6L 

or  H  =  l.l/L)  partly  to  correct  for  dynamic  effects  and  partly 

because  it  is  thought  than  an  L/20  wave  is  too  severe  for  the 

very  large  ships.   Notwithstanding,  this  empirically  derived 

procedure  was  considered  entirely  satisfactory  until  relatively 

recently  and  is  still  in  prevalent  use  inasmuch  as  it  gives 

conservative  estimates. 

The  appearance  of  large  ocean-going  tankers  and  bulk 
carriers  and  novel  ship  types  gave  impetus  towards  the  search 
for  a  "rational"  procedure  for  ship  structural  design.   It 
was  realized  that  the  application  of  the  standard  bending 
moment  calculations  for  these  ships  lead  to  unrealistic  re- 
sults.  At  the  same  time  new  knowledge  relating  to  ship 
structural  design  had  accumulated  and  experimental  and  comp- 
utational facilities  to  advance  this  knowledge  have  improved 
tremendously . 


Many  marine-oriented  societies  and  organizations 
undertook  basic  research  to  improve  the  accuracy  of  pre- 
dictions of  sea  loads  and  motions,  and  the  development  of  a 
rational  design  method  applicable  for  all  ship  types.   Full 
scale  and  model  tests  data  were  collected.   The  random 
nature  of  both  the  loads  on  a  ship  and  its  capacity  to  re- 
sist these  loads  was  recognized,  so  that  probability  theory 
was  utilized  to  develop  statistical  models  for  the  design  of 
ships  [1-14]*.   The  advantages  of  a  probabilistic  approach 
over  the  traditional  deterministic  calculations  are  discussed 
in  reference  [7]. 

In  the  conventional  method,  the  primary  strength  of 
a  ship  is  based  solely  on  the  longitudinal  stresses  caused 
by  vertical  bending,  or  bending  about  the  horizontal  neutral 
axis.   The  effects  of  lateral  bending,  or  bending  about  the 
vertical  neutral  axis,  is  neglected.   However,  for  ships 
travelling  oblique  to  the  waves,  the  lateral  loads  reach  a 
magnitude  equal  to  or  greater  than  the  vertical  bending 
moment.   The  combined  effect  of  lateral  and  longitudinal 
bending  may  be  critical  for  some  ship  types. 

A  requisite  for  an  integrated  probabilistic  design 
procedure  is  a  more  precise  knowledge  of  the  loads  acting  on 


*The  numbers  in  brackets  refer  co  references  listed 
in  the  Appendix 
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a  ship.   This  paper  represents  a  small  portion  of  the  total 
effort  to  completely  define  the  structural  loads  experienced 
by  a  ship  in  a  seaway.   It  is  an  investigation  of  the  magnitude 
of  combined  horizontal  and  vertical  bending  in  irregular  seas. 
Due  to  their  relative  phase  difference,  the  amplitudes  of  the 
vertical  and  horizontal  moments  generally  are  not  acting  at 
the  same  time . 

In  the  analysis,  the  loads  are  treated  as  random 
variables.   Consideration  of  the  total  longitudinal  stress  at 
the  deck  leads  to  the  definition  of  an  equivalent  combined 
bending  moment.   Two  methods  are  developed  to  generate  the 
statistical  parameters  of  the  combined  load  so  that  the  usual 
statistical  methods  to  predict  long-term  extreme  values  of  the 
load  can  be  applied.   An  example  application  is  then  made  for 
a  large  tanker  which  had  been  previously  instrumented  so  that 
full  scale  data  for  the  vertical  bending  stresses  experienced 
in  service  are  available  for  purposes  of  comparison.   Such  a 
comparison  serves  to  check  yet  another  time  the  validity  of 
computational  techniques  to  simulate  ship  motion  and  loads  in 
waves,  and  in  particular  the  MIT  5-D  seakeeping  program.   The 
predicted  moment  responses  are  compared  and  analyzed  for  trends , 


Chapter  II 

PROCEDURE 

II. 1   Determination  of  the  Equivalent  Combined  Wave  Bending 
Moment 

A  ship  travelling  at  an  angle  to  the  waves  will  be 

subjected  to  unsymmetrical  bending  such  that  the  stress  at 

the  deck  is  not  uniform  but  assumes  a  maximum  value  at  one 

edge,  as  shown  in  figure  II. 1.   This  primary  stress  may  be 

thought  of  as  the  combined  effect  of  vertical  and  horizontal 

bending  and,  to  a  much  lesser  extent,  torsion.   Previous 

studies  [15,  16]  have  established  that  torsional  moments  are 

relatively  small,  never  exceeding  10%  of  the  vertical  moment 

in  any  given  condition.   In  this  study  the  effects  due  to 

torsion  are  considered  negligible,  but  it  is  recognized  that 

torsion  loading  may  be  critical  for  some  ships,  as  those 

with  large  hatch  openings  [17].   The  contribution  to  the 

ship  primary  stress  may  be  included,  if  so  desired,  in  a 

manner  analogous  to  the  method  described  here  [18]. 

The  stress  at  the  deck  edge  is  given  by 


°T  =  °V  +  °H 


+  c 


+ 


Vertical  Bending  Stress 
(sagging) 


Horizontal  Bending  Stress 


Combined  Stress 

Figure  II .  1 

The  Longitudinal  Stress  Distribution 
on  a  Ship  Hull  Girder 


where 


o      =  combined,  or  total,  stress 


o     =  stress  due  to  vertical  bending 


o„  =   stress  due  to  horizontal  bending 


In  terms  of  moments,  one  can  write 


=  *V  +  MH 
T    *V    ^H 


where 


ML  =  vertical  bending  moment 


M„  =  horizontal  bending  moment 
H 


ivv 
k T  =  JJ  =  vertical  section  modulus 
V    z 


k   =   zz  =  horizontal  section  modulus. 
h    y 

multiplying  both  sides  of  the  equation  by  k^,  one  gets 
°TkV  =  "V  +  ^  MH 


If  one  defines  an  equivalent  combined  bending  moment 

^  =  0rrK7,    and  the  section  modulus  ratio  K  =  ■*— -,    then  one 

gets  the  relation 


MT  =  My   +  KMR  (1) 

It  is  emphasized  here  that  one  cannot  simply  add 
the  maximum  values  of  the  vertical  and  horizontal  moments  to 
get  the  maximum  combined  effect.   In  general  the  moments  at 
any  point  along  the  length  of  the  ship  are  not  acting  in 
phase.   Furthermore,  it  may  be  of  interest  to  specify  the 
complete  distribution  of  the  primary  stress  at  the  deck. 
Within  the  limits  of  the  beam  theory,  this  is  done  by  getting 
the  stress  at  one  other  point  across  the  ship's  deck.   A 
convenient  reference  point  is  at  the  ship's  center line  where 
the  stress  due  to  lateral  bending  is  zero  for  most  ships  due 
to  port-and-starboard  symmetry,  i.e.,  a     =   av  at  the  centerline. 

In  the  statistical  approach  to  sea  loads,  the  ex- 
citing force  is  assumed  to  be  a  stationary,  narrow-band  and 
gaussian  process  with  zero  mean  over  a  short  period  of  time 
(e.g.  30  minutes).   The  ship  is  considered  to  be  a  linear 
system  so  that  the  response  is  also  stationary,  narrow-band 
and  gaussian  with  zero  mean  over  the  same  period  of  time. 
Then  the  amplitudes  follow  the  Rayleigh  distribution  [19,  20]. 
Over  a  long  period  of  time,  the  assumption  of  a  stationary 
excitation  process  cannot  be  made.   However,  it  was  found 
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that  the  amplitudes  follow  closely  the  Weibull  distribution 
[8,  9,  10].   It  is  shown  that  the  Rayleigh  is  a  special  case 
of  the  Weibull  distribution  given  by: 


fx(x)  =  (s/k)(x/k)S  1   exp-(x/k)S     X10  (2) 


Fx(x)  =  1  -  exp-(x/k)S        X>0       (3) 


wher  fy(x)  and  F  (x)  are  the  probability  density  and  distrib- 
ution functions  respectively,  and  s  and  k  are  parameters.   To 
get  the  Rayleigh  distribution,  one  sets  s=2  and  k  =  /E  where 
E  is  the  mean  square  value  of  the  amplitude  of  the  response. 
Furthermore  it  was  found  that  the  best-fit  Weibull  curve  is 
approximately  the  exponential  distribution  obtained  from 
equations  (2)  and  (3)  by  setting  s  =  1  and  k  =  m  where  m  is 
the  mean  value  of  the  response  amplitude.   The  statistics  of 
the  ship  bending  moment  response,  therefore,  can  be  specified 
if  the  mean  square  value,  or  variance,  of  the  process  is 
known . 

The  usual  method  to  obtain  the  variance  of  the  re- 
sponse in  random  seas,  attributed  to  St.  Denis  and  Pierson, 
consists  of  specifying  the  energy  density  spectra  of  the 
ocean  waves  and  obtaining  the  ship  response  amplitude  opera- 
tor (RAO) ,  or  the  response  per  unit  amplitude  of  regular 
wave,  by  model  tests  or  some  computational  procedure  using 
the  strip  theory  approach.   From  the  theory  of  linear  systems, 
the  spectral  density  of  the  response  c£'n  be  found  from  the 


relation  [21] 


SR(u>)  =  |H(u>)  |  2S(o)) 


where  S  (to)  =  spectral  density  of  the  response 


H(co)  =  ship  response  amplitude  operator  (RAO) 

S(w)  =  spectral  density  of  the  seaway. 

The  area  under  the  response  spectrum  is  the  variance 
of  the  process  and  represents  the  sum  of  the  responses  to  the 
infinitely  many  sine  waves  that  are  thought  to  compose  the 
actual  sea. 

Two  methods  used  to  obtain  the  variance  of  the  wave 
bending  moments  by  computer  simulation  techniques  are  ex- 
plained below. 

1)   First  Method 

Equation  (1)  can  be  written  showing  the  moments  as 
a  continuous  random  function  of  time,  t: 


M 


T(t)  =  MyU)  +  KMH(t) 
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Using  the  assumption  that  the  sea  can  be  represented  by  in- 
finitely many  harmonic  waves  of  small  amplitudes,  continuous 
frequency  and  random  phase,  one  can  write  the  vertical  and 
horizontal  moment  responses  as: 


M,7(t)  =  V  cos  O  t  +  £  )  ,  and 
v  n    v  n     n; 


M  (t)  =  H  cos  (u  t  +  6  ) 
H       n     n     n 


where 


V  =  amplitude  of  the  vertical  moment  response  to  the 
nth  wave 


H  =  amplitude  of  the  horizontal  moment  response  to 


nth  wave 


a)  =  frequency  of  the  nth  wave 

e   =  phase  angle  of  the  vertical  moment  with  respect 


to  a  wave  with  its  crest  amidship 


9   =  phase  angle  of  the  horizontal  moment  with  respect 
n 


then 


M 


to  a  wave  with  its  crest  amidship 


m(t)  =   A-  V  cos  (to  t  +  £  )  +  K  Lh  cos(ft)  t  +  6  ) 
T      n-1   n    v  n     n'     n=l  n    v  n     n' 


or 


on 

m(t)  =   A  [v  cos(ft)  t  +  e  )  +  KH  cos(u  t  +  6  ) 
T      n=lL  n     n     n      n    ^  n     n'J 


M 


which  can  be  rewritten  as 


11 


M„(t)  =  I 
n=l 


(v  cos  e  +  KH  cos  6  )  cos  u  t 
x  n     n     n     n'      n 


Let 


Then, 


-  (V  sin  e   +  KH  sin  9  )    sin  u  t 
v  n     n     n     n '  n  i 


X   =  V  cos  e  +  KH  cos  0 
n    n     n     n     n 


Y  =  V  sin  e  +  KH  sin  6 
n    n     n     n     n 


n=l  L 


X  cos  to  t  +  Y  sin  w  t 
n     n     n     n 


(4) 


Since  M  (t)  is  also  a  continuous  random  function  of  time,  it 
can  also  be  represented  by: 


Mm(t)  =   Y  Z  cos  (a)  t  +  e  ) 
T        i  n    v  n     n; 
n=l 


(5) 


A  comparison  of  equations  (4)  and  (5)  shows  that 


-1  Y  ?     2     2 

6   =  tan     n,  and   Z   =  X   +  Y 

n         ~v~  n     n    n 

n 
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Then 


Z         =  V        +   K   H        +    2KH  V   cos(e      -    6    )  (6) 

n  n  n  nnvnrr 


Equation  (6)  gives  the  amplitude  of  the  combined 
moment  due  to  vertical  and  horizontal  bending,  with  the 
phase  difference  in  the  occurence  of  the  component  peaks 
taken  into  consideration.   The  MIT  5-D  seakeeping  program 
calculates  the  ship  motions  and  the  vertical  and  horizontal 
moment  responses  to  regular  waves  (RAO) .   A  slight  modifica- 
tion of  the  program  to  obtain  the  equivalent  combined  moment 
response  operator  was  made  by  use  of  equation  (6).   The  RAO's 
were  then  used  in  the  usual  power  spectral  analysis  in  the 
frequency  domain  to  obtain  the  variance  of  the  responses  in 
irregular  seas.   Comparison  of  predicted  and  experimental  re- 
sults have  shown  generally  good  agreement  [15,  16,  21,  22,  23] 
and  have  led  to  widespread  acceptance  of  the  superposition 
method  of  St.  Denis  and  Pierson. 

2)   Second  Method 

Considering  the  wave  loads  as  random  variables,  it  is 
generally  accepted  that  the  bending  moments  M^  and  M   are  sta- 
tionary, narrow-band,  gaussian  processes  with  zero  means  for 
short  periods  of  time.   Then  M   is  also  stationary,  narrow-band 
and  gaussian  with  zero  mean,  i.e.,  in  statistical  notations: 

EM    =  E  ^  +  KE  M    =0 
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Furthermore, 


E 


M„ 


=  E 


K  +  mnY 
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=    E 

H] 

+    K2E 

M2" 
H 

+    2KE 

^ 


M 


H 


but  the  covariance  is  defined  as 


E 


MyM 


H 


A/hHK"  ^)2]e[(mh  "  ™uf}}h 


where  m^  =  m  =  mean  values  =  0 


fL„  -   correlation  coefficient  of  the  vertical  and 
horizontal  bending  moments 


therefore, 


MI 


A 


7 

+    K   E 


R] 


+ 


2k/3vh(eK 


M 


H 


(7) 


For  a  zero-mean  process,  the  mean  square  is  equal  to  the 
variance. 
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As  in  the  first  method,  the  value  of  the  section 
modulus  ratio,  K,  must  be  known  for  a  given  ship  or  esti- 
mated during  preliminary  design.   An  attempt  was  made  to 
examine  the  section  modulus  ratio  of  several  sizes  of  tankers 
to  determine  if  there  was  a  trend  with  the  growth  in  the  size 
of  these  ships.   As  shown  in  figure  II. 2,  the  points  are 
widely  scattered  with  no  apparent  trend.   Since  the  value  of 
K  depends  on  the  construction  of  the  ship  (framing  system, 
scantlings,  etc),  perhaps  the  scatter  of  points  is  indica- 
tive of  the  present  uncertainty  about  both  sea  loads  and  the 
strength  of  ships,  and  the  rather  arbitrary  way  a  "safety 
factor"  is  arrived  at. 

Likewise,  the  value  of  the  correlation  coefficient, 
/>  ,  of  the  vertical  and  horizontal  moments  must  be  deter- 
mined.  Implicit  in  the  use  of  equation  (7)  is  the  assumption 
that  an  average  value  of  /O   already  takes  into  account  the 
phase  difference  of  the  peaks  of  M^  and  M   in  all  conditions 
at  sea. 

With  this  method,  the  bending  moment  response  ampli- 
tude operators  may  be  determined  by  experiment  or  computer 
techniques.   Spectral  analysis  could  then  be  performed  to 
obtain  the  variances  of  M^  and  M  ,  and  finally  the  combined 
moment  M 
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II. 2   Long-Term  Wave  Load  Analysis 

After  the  statistics  of  the  short-term  response  are 
determined  for  the  range  of  sea  states  that  make  up  the 
actual  seaway,  the  expected  long-term  loads  may  be  predicted. 
Several  methods  to  obtain  the  long-term  distribution  of  wave 
loads  and  the  prediction  of  extreme  values  have  been  proposed 
[6,  9,  11,  18,  23].   Mansour  [18]  has  described  procedures 
for  the  analysis  of  failure  under  extreme  loads  under  two 
conditions: 

1.  Short-term  analysis  -  when  the  route  of  the  ship 
is  not  known  and  a  "severest  condition"  criterion 
is  imposed. 

2.  Long-term  analysis  -  If  the  ship's  route  is  known 
and  if  that  route  is  more  or  less  permanent. 

In  this  study,  a  procedure  is  discussed  to  determine 
the  value  of  the  parameter  k  =  m  of  equations  (2)  and  (3) 
for  a  ship  travelling  oblique  to  waves  in  all  sea  conditions 
described  by  the  significant  wave  height  and  period  (h,T)  as 
in  [24].   The  following  assumptions  are  made: 

1.   wave  headings  are  uniformly  distributed  over  all 
directions . 


17 


2.   The  speed  is  essentially  determined  by  the  state 
of  the  sea.   The  speed-significant  wave  height  relation  is 
adopted  from  [14]  as  follows: 


Speed,  V   knots 


significant  wave  height,  h   ft, 


16 

12 

6 


less  than  12.5 
12.5  to  21.5 
greater  than  21.5 


The  procedure  followed  is  described  below: 

1.  Obtain  the  RMS  response  to  all  sea  states  (i.e., 
range  of  sea  spectra)   in  all  conditions  of  speed 
and  wave  headings  by  some  analytical  method  as 
the  "MIT  5-D  seakeeping  program." 

2.  For  each  speed,  get  the  average  of  the  RMS 
response  over  all  wave  directions. 

3.  Use  Bretschneider ' s  two  -parameter  family  of  sea 

spectra  to  get  the  response  to  any  sea  described 

by  a  particular  (h,T)  pair.   This  is  most  easily 

done  by  scaling  the  responses  to  the  sea  spectra 

of  step  (1) .   We  call  this  matrix  of  the  response 

R  with  elements  r. .. 

ID 


4.   For  each  wave  height  in  [24],  determine  the 

probability  (i.e.,  relative  frequency  of  occurence) 
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of  each  observed  period.   We  call  the  resulting 
matrix  of  probabilities  P  with  elements  p. . . 


5.   Determine  the  average  RMS  response  for  each 
significant  height  by  the  relation 


m.  =  7  r .  .  x  p.  . 
i    j   ID     ID 


where  the  subscripts  i  and  j  refer  to  rows  and 
columns  of  a  matrix.   The  resulting  column 
matrix  is  called  M. 

6.  Determine  the  probability  of  occurrence  of  each 
significant  wave  height.   We  call  this  column 
matrix  Q  with  elements  q. . 

7.  The  mean  or  expected  value  of  the  RMS  of  the 
response  for  a  particular  area  is  given  by 


mean  RMS  =  J   m.  x  q. 

l 


Then  the  desired  parameter  m  =  1.25  x  (mean  RMS 
of  the  response)  in  the  case  of  a  short-term 
analysis  [18].   For  the  long-term  analysis,  a 
time  factor  is  additionally  introduced  to  get  the 
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weighted  mean  of  the  RMS  of  the  response  for  all 
areas  in  the  ship's  route. 

Steps  4-6  is  equivalent  to  getting  the  probability 
of  occurrence  of  each  (h,T)  pair  in  any  one  area  being  de- 
scribed.  If  we  call  this  matrix  P'  with  elements  p! . ,  then 

ID 


mean  RMS  =  7  Y  r .  .  x  p'.  . 
h  h  ii  ri] 
i  D    J      J 


The  data  in  [24]  have  been  computerized  and  could 
be  used  in  the  above  calculations. 
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II. 3   Details  of  the  Procedure 

The  methods  for  wave  load  analysis  and  prediction,  as 
explained  in  the  foregoing  sections  of  this  chapter,  were 
applied  in  the  analysis  of  a  large  tanker.   The  UNIVERSE 
IRELAND  was  one  of  five  vessels  of  different  sizes  instrumented 
for  full-scale  measurement  of  stresses  experienced  in  actual 
service.   However,  the  strain  gage  arrangement  on  deck  was 
such  that  the  stress  due  to  lateral  bending  was  cancelled  out, 
so  that  only  vertical  bending  stresses  were  recorded.   There- 
fore, the  results  derived  for  the  response  of  primary  interest, 
the  combined  moment,  could  not  be  verified.   The  description 
and  analysis  of  the  results  of  the  3^-year  project  sponsored 
by  the  American  Bureau  of  Shipping  (ABS)  is  reported  in  [25] 
and  provide  the  basis  for  comparison  with  the  analytically 
derived  values  of  vertical  bending  stress  at  amidship. 

The  general  characteristics  of  the  UNIVERSE  IRELAND 
are  listed  below: 


Full  load  displacement,  tons 375,811. 

Length  overall,  ft  1135.17 

Length  between  perpendiculars,  ft   ...      1083. 

Breadth,  ft   174.87 

Depth,  ft 10  5. 

Design  draft  (keel)  ,  ft 81-0-5 

Block  coefficient  (kwl) 0.86 

Section  modulus,  top,  in2-ft  566,794. 

Section  modulus,  side,  in2-ft   617,350. 


The  ship's  lines  and  other  data  pertaining  to  the 
ship  necessary  to  generate  the  input  for  the  computer  analysis 
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were  obtained  from  the  ABS .   It  was  necessary  to  use  an 
approximate  lightship  weight  distribution  because  of  an 
apparent  discrepancy  in  the  distribution  curve  provided. 
Since  the  still  water  bending  moment,  at  least,  is  very 
sensitive  to  the  distribution  of  the  weight,  a  probable 
error  was  anticipated  in  the  calculated  values.   For  ex- 
ample, the  maximum  still  water  bending  moment  calculated 
by  the  computer  exceeded  the  value  determined  by  the  ABS 
by  about  11.5%. 

The  MIT  5-D  seakeeping  program  is  fully  described 
in  [26].   Only  minor  modifications  were  necessary  to  do  the 
quasi-static  combined  bending  moment  analysis.   The  re- 
sponse to  long-crested  irregular  seas  were  used  in  conjunc- 
tion with  an  auxiliary  program  to  obtain  the  responses  to 
short-crested  seas.   This  program  uses  the  cosine-squared 
spreading  function 


s(/<)  =\   cos2(/*) 


where/tt=  angle  between  the  wind  direction  and  the  wave 
components.   The  principal  wind  directions  were  given  ex- 
actly the  same  values  as  the  wave-to-course  headings. 

The  computer  was  run  for  ship's  speed  of  6,  12  and 
16  knots,  nine  wave  headings  from  0   to  180   in  increments  of 
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22.5  ,  and  nine  sea  states  described  by  the  significant 
height  and  period.   The  computer  output  consisted  of  the 
root  mean  square  of  the  response,  the  significant  height 
(H/?)  and  the  average  of  the  l/10th  highest  (H  no)      of  the 
response  measured  peak-to-peak  for  all  combinations  of  the 
three  speeds,  nine  wave  directions  and  nine  sea  states. 
Since  it  was  desired  to  compare  the  values  of  the  vertical, 
horizontal  and  combined  wave  moments,  output  was  obtained 
for  all  three  responses.   Furthermore,  to  be  as  accurate 
as  possible,  the  three  responses  were  calculated  for  twenty 
stations  specified  along  the  length  of  the  ship,  even  if 
eventually  only  the  responses  amidships  were  utilized  in 
the  analysis. 

The  sea  data  were  carefully  chosen  to  match  the 
actual  seaway  experienced  by  the  ship  during  the  instru- 
mentation period.   This  was  done  by  obtaining  sea  data  from 
[24]  for  the  areas  lying  along  the  service  route  of  the 
UNIVERSE  IRELAND,  which  made  regular  runs  between  the  Persian 
Gulf  and  Western  Europe  by  way  of  the  cape  of  Good  Hope  [25]. 
Although  the  seakeeping  program  can  accept  and  use  an  actual 
sea  spectrum  or  simulate  developing  or  decaying  seas  by  the 
Bretschneider  two-parameter  spectrum,  the  Pierson-Moskowitz 
one-parameter  family  for  fully  developed  seas  was  chosen. 
Nine  sea  states  with  significant  heights  ranging  from  3.6  feet 
to  43.3  feet  (corresponding  to  observed  periods  of  1.5  seconds 
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to  18.5  seconds)  were  used.   Spectral  ordinates  for  forty 
frequencies  were  calculated  while  twenty  five  points  were 
used  to  describe  the  response  amplitude  operators  covering 
the  wave  frequency  range  0.1  <  -  <  9.0  with  most  of  the 

Li  — 

points  spaced  closely  around  the  expected  spectral  peak  at 
-  =  1.0.   A  is  the  wavelength  and  L  is  the  length  of  the 

Li 

ship.   The  response  to  regular  waves  (RAO's)  were  not 
printed  in  order  to  save  on  the  cost  of  running  the  computer. 


Chapter  III 
RESULTS 
The  following  general  areas  of  study  were  explored: 

1.  Comparison  of  the  derived  results  with  those 
in  other  published  works  on  the  subject. 

2.  Comparison  of  the  variance  of  the  combined 
moment  obtained  from  the  two  methods  described 
in  the  previous  chapter. 

3.  Comparison  of  theoretical  calculations  for  the 
vertical  bending  stress  with  full-scale 

measurements. 

4.  Examination  of  bending  moments  trends  and  wave 
load  prediction. 

The  results  of  the  above  studies  are  presented  in 
the  following  sections. 
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III.l   Comparison  of  Results  with  Other  Published  Works 

The  calculated  value  of  the  root-mean- square  (RMS) 
of  the  bending  moment  responses  are  plotted  in  figures  III.l 
to  III. 4  as  a  function  of  the  wave  headings.   For  contrast 
this  was  done  for  the  responses  to  long-crested  and  short- 
crested  fully  developed  seas  with  significant  wave  heights 
of  7.4  ft.,  19.2  ft.  and  4  3.3  ft.   It  was  intended  to  show 
the  response  to  a  low,  moderate  and  high  sea  states.   Since 
only  the  value  of  the  horizontal  moment  in  "realistic"  seas 
are  of  interest,  the  RMS  of  the  response  were  plotted  only 
for  short-crested  seas.   The  graphs  were  prepared  for  a 
representative  ship's  speed.   The  general  behavior  of  the 
response  did  not  change  with  speed. 

Faltinsen  [22]  describes  the  comparison  between 
theory  and  experiments  of  tne  wave-induced  horizontal  and 
vertical  forces  and  moments  of  a  series  60  hull  form  with  a 
block  coefficient  of  0.80  and  a  Froude  member  of  0.15  moving 
in  regular  waves  with  different  headings  at  an  angle  to  the 
direction  of  the  waves.   The  predictions  made  by  the  new 
strip  theory  developed  by  Salvesen,  Tuck,  and  Faltinsen, 
which  is  used  in  the  MIT  5-D  seakeeping  program  proved  very 
satisfactory.   The  theory  is  described  by  Salvesen  et  al,  in 
[27].   For  vertical  bending  moments,  the  theory  predicted 
higher  values  than  the  experiments  for  -  <  0.7  in  quartering 
and  following  waves  and  lower  values  for  X/L  <  0.6  in  bow 
waves.   This  trend  may  also  be  observed  by  comparing  the  RMS 
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values  of  ML^  in  short-crested  seas  for  wave  headings  of  0° 
(following  seas)  and  180   (head  seas)  in  figures  III. 2,  3 
and  4.   For  horizontal  bending  moments,  appreciable  diff- 
erence existed  between  theory  and  experiments  only  for  beam 
seas. 

One  of  the  early  investigations  of  the  effect  of 
horizontal  moments  in  oblique  seas  is  described  by  Numata 
[16].   It  was  established  in  this  study  that  lateral  moments 
reach  values  that  are  substantial  and  that  peak  bending 
moment  exists  in  the  region  of  the  speed  for  synchronous 
pitching  and  heaving  motions.   Wahab  [15]  shows  that  the 
importance  of  the  horizontal  bending  moment  increases  with 
increasing  ship  length.   For  the  sea  condition  used  in  his 
investigation,  the  significant  midship  horizontal  moment  for 
a  ship  of  200  meter-length  was  46  per  cent  of  the  significant 
vertical  moment  while  for  a  400  meter-length  ship  the  per- 
centage increased  to  69.   The  concept  of  effective  wavelength, 
L  ,  described  by  Numata  was  helpful  to  explain  the  results 
obtained  from  the  computer  analysis.   However,  only  one  sea 
spectrum  was  used  in  his  study  and  some  conclusions  were 
reached  that  are  not  necessarily  true  for  all  conditions  at 
sea  for  any  given  ship.   Inasmuch  as  the  RAO  of  the  UNIVERSE 
IRELAND  was  not  recorded  by  the  computer,  the  figures  in 
[15]  proved  very  useful.   The  ship  studied  in  that  work  is 
the  same  as  the  model  described  by  Faltinsen  and  has  general 
characteristics  similar  to  the  UNIVERSE  IRELAND.   The  respcnse 
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of  the  model  to  regular  waves  could  be  assumed  to  be  sub- 
stantially identical  in  shape  and  trends.   The  calculation 
of  amidship  forces  and  bending  moments  in  [15]  also  made 
use  of  a  single  spectrum  representing  a  severe  sea. 
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III. 2   Comparison  of  the  Calculated  Values  of  the  Combined 
Moment  Response 

Table  III-l  is  a  summary  of  the  results  obtained 
for  the  RMS  value  of  the  combined  moment  response  for  the 
vessels'  service  speeds  of  16  and  6  knots.   The  first  set 
of  values  are  those  calculated  from  the  RAO  developed  by  the 
computer  by  use  of  equation  (6)  of  the  first  method.   The 
second  set  of  values  were  calculated  from  equation  (7)  of 
the  second  method  using  the  value/3   =  0.32  obtained  from 
[18]  as  reprinted  from  the  published  result  of  a  Det  Norske 
Veritas  Research  Report.   The  reprinted  source  is  shown  as 
figure  III. 5.   The  values  of  P        represent  the  correlation 
coefficient  of  the  most  probable  largest  values  of  the 
vertical  and  horizontal  moment  for  one  ship  life  of  about 
twenty  years  in  the  North  Atlantic.   They  were  meant  to  be 
reference  values  and  not  for  use  in  design.   The  difference 
in  M  values  is  shown  for  each  sea  state.   It  is  immediately 
apparent  that  the  results  obtained  by  both  methods  are  not 
very  different.   The  third  set  of  values  was  calculated  also 
by  use  of  equation  (7)  but  using  a  value  of  /^   obtained  by 
averaging  the  responses  in  short-crested  seas  by  the  first 
method  for  all  wind-wave  directions  and  for  three  represent- 
ative spectra  for  low,  moderate  and  high  seas.   The  mean 

value  of  p.7ZJ   obtained  in  this  manner  was  0.532,  significantly 
'  VH 

higher  than  the  ISSC  value.   The  difference  in  the  results  for 
M   is  also  indicated  for  comparison. 
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Figure    III. 5 

Correlation    Coefficient    Between    Vertical    Bending    Moment 
and    Horizontal    Bending    Moment 
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III. 3   Comparison  of  the  Calculated  Vertical  Bending  Stress 
and  Full-scale  measurements; 

Strain  gages  were  installed  on  the  port  and  starboard 
sides  at  amidship  and  several  other  stations  along  the  length 
of  the  UNIVERSE  IRELAND  to  measure  stresses  on  the  deck.   The 
arrangement  was  such  that  only  vertical  bending  stresses  were 
measured.   Therefore,  though  the  combined  moment  could  not  be 
verified,  a  comparison  between  full  scale  data  and  the  com- 
puter-derived values  was  possible.   This  comparison  by 
weather  group  as  specified  by  wind  measurements  (Beaufort 
scale)  is  shown  in  figure  III. 6.   A  summary  is  also  pre- 
sented in  table  III. 2.   The  conversion  from  bending  moment 
to  stress  is  done  by  the  relation:   Bending  Moment  =  Stress  x 
Section  Modulus.   The  values  compared  are  the  average  RMS 
values  of  the  peak-to-peak  stress  for  each  wind  scale.   The 
full  scale  data  were  read  off  figure  11  of  [25].   The  con- 
version from  wind  to  wave  was  done  by  the  graphical  relation 
shown  in  figure  III. 7  [23].   This  conversion  is  summarized 
in  table  III. 3. 

The  use  of  nine  sea  spectra  with  each  significant 
wave  height  enables  the  calculation  of  the  standard  devia- 
tion of  the  predicted  values  of  the  vertical  stress.   For 
comparison  with  the  measured  data,  these  standard  deviations 
had  to  be  modified  to  obtain  values  based  on  wind  measure- 
ments.  This  was  accomplished  by  the  relation  [23]: 


2      2 

Sl   "  V  + 


S  2  -  AH 
J      2 
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2 

tan  6 


where    S.  =  standard  deviation  of  the  response  with  respect 
to  the  wind 


S„  =  standard  deviation  of  the  response  with  respect 
to  waves 


S_  =  standard  deviation  of  the  wave  height  with  re- 
spect to  the  wind,  as  shown  in  figure  III. 8 
from  C23 J. 


AH  =  wave  height  increment 

8  =  average  slope  of  the  response  vs.  significant 
height  plotted  in  figure  III. 9. 

For  comparison,  table  III. 4  of  the  standard  deviations  fo: 
all  the  cases  studied  is  drawn  up. 
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III. 4   Examination  of  Bending  Moment  Trends  and  Wave  Load 
Prediction 


Since  the  traditional  primary  strength  calculations  use 
the  vertical  bending  moment  in  a  long-crested  sea  of  arbirtary 
height  as  a  basis  for  ship  load,  it  is  of  interest  to  compare 
the  values  of  the  vertical  bending  response  in  long-crested  head 
seas  with  the  combined  moment  in  short-crested  seas  for  the 
higher  range  of  the  wave  significant  heights.   Figures  III. 10 
and  III. 11  are  drawn  for  sea  spectra  representing  fully  developed 
moderate  and  severe  seas,  respectively  moderate  seas,  the  com- 
bined moment  for  head  seas  exceeded  the  long-crested  vertical 
moment  by  about  thirty  sever  per  cent  (37%)  and  in  severe  seas, 
the  difference  was  sixteen  per  cent  (16%) . 

Similarly,  it  was  deemed  desirable  to  compare  the  rel- 
ative values  of  the  vertical  and  horizontal  moments  for  all 
conditions  at  sea.   This  is  done  in  figures  III  2,  4  and  6  for 
the  responses  in  short-crested  seas.   One  can  see  that  the 
horizontal  moment  is  not  small  compared  to  vertical  moment; 
and  in  low  to  moderate  seas,  the  former  is  at  least  as  great  as 
the  latter.   The  above  mentioned  comparison  of  the  response  in 
long-crested  and  short-crested  head  seas  reflect  the  relative 
contributions  of  the  vertical  and  horizontal  moments  to  the 
total  load. 

Table  III. 5  illustrates  the  change  with  speed  in  the 
bending  moment  response  for  low,  moderate  and  high  seas.   The 
tabulated  values  are  the  mean  RMS  peak-to-mean  responses  in 
short-crested  seas  averaged  over  all  wind-wave  directions. 
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In  order  to  illustrate  the  procedure  to  obtain  the 
parameter  k  for  the  long-term  distribution  of  the  wave  re- 
sponse, a  "short-term"  analysis  was  done  for  the  worst 
condition  in  the  ship's  route  (area  7  in  [24]).   This  is  the 
North  Atlantic  in  the  months  of  December  to  February  and  the 
tabulated  data  of  [24]  were  used.   The  results  obtained  are 
given  as  tables  III. 6  to  III. 9  for  the  combined  moment.   For 
this  analysis,  the  result  is 

mean  RMS  =  367820  ft-tons 

Since  this  is  a  "short-term"  analysis,  one  can  then  say 
m  =  1.25  (mean  RMS)  =  459775  ft. tons. 

The  equivalent  values  for  the  vertical  and  horizon- 
tal moments  were  also  computed.   These  are: 

mean  RMS  response,  M^  =  19  730  7  ft-tons 


mean  RMS  response,  M„  =  2  310  4  3  ft-tons 


Note  that  the  expected  value  for  the  horizontal  moment  is 
larger  because  the  ship  encounters  more  low-to-moderate 
seas  in  any  length  of  time  than  severe  seas. 
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The  mean  value  of  the  RMS  of  the  response  was  also 
calculated  for  all  the  areas  in  the  route  of  the  ship,  i.e., 
a  "long-term"  analysis  was  conducted.   Both  IVL^  and  M  were 
calculated  and  the  results  are  summarized  in  table  III. 10. 
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Table  III. 6 

Mean  RMS  of  the  Combined  Moment  Response  in  Short-Crested  Seas 
Averaged  Over  Uniformly  Distributed  Wind-Wave  Directions 


Velocity 

Significant 

RMS  Mv 

RMS  MH 

RMS  MT 

(knots) 

Wave  Height 
(ft) 

(ft-tons) 

(ft-tons) 

(ft-tons) 

3.6 

13386.3 

21411.2 

26762.9 

7.4 

54212.8 

79641.4 

115467.5 

10.7 

115913.9 

159110.8 

240620.1 

14.6 

212240.4 

253534.2 

405516.0 

V=16 

19.2 

334744.1 

346676.0 

586564.9 

24.3 

463645.8 

426305.6 

756493.1 

30.0 

591299.8 

488621.4 

911865.7 

36.4 

712415.3 

537815.8 

1050905.6 

43.3 

924782.1 

573471.8 

1,167,628.3 

3.6 

11689.7 

19464.6 

24197.6 

7.4 

51434.1 

78228.2 

112,590.6 

10.7 

112157.7 

160,233.2 

238895.4 

14.6 

206043.0 

258953.5 

403795.1 

V=12 

19.2 

327646.2 

358380.9 

585867.4 

24.3 

456173.6 

442364.9 

757601.8 

30.0 

583590.6 

509395.5 

914721.5 

36.4 

705103.0 

561731.7 

1055359.7 

43.3 

810988.6 

597580.8 

1172221.1 

3.6 

9831.6 

17268.6 

22238.0 

7.4 

49947.6 

77878.1 

112551.8 

10.7 

110104.1 

157927.6 

236043.6 

14.6 

201175.2 

258401.3 

398461.4 

V=6 

19.2 

320621.7 

360908.4 

580125.1 

24.3 

449758.8 

449322.7 

756008.3 

30.0 

696331.0 

518255.7 

913820.7 

36.4 

696331.0 

573820.7 

1056466.6 

43.3 

801856.7 

6] 5268. 6 
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Table  III. 10 

Mean  RMS  Value  of  the  Response  for  all  Areas  in  the 
Route  of  the  UNIVERSE  IRELAND 


Area 

E[RMS]of  M^  Response 

E[RMS]of  Mp  Response 

Time  Factor 

7 

154254 

299294 

2 

11 

126213 

246724 

2 

18 

109398 

218546 

1 

20 

66961 

138544 

1 

23 

110971 

220161 

1 

28 

80454 

170798 

1 

29 

92992 

189706 

1 

34 

108754 

203516 

1 

41 

146576 

282585 

2 

42 

123456 

246480 

2 

45 

161230 

260578 

2 

Time  Weighted  Average: 
(RMS)V  =  124,562  ft-tons 


(RMS)   =  517,450  ft-tons 


m^  =  155,702  ft-tons 
iiL  =  646,813  ft-tons 


—  =  4.15 
*V 


Chapter  IV 

DISCUSSION  OF  RESULTS 

IV. 1   Evaluation  of  the  analytical  Results 

The  results  obtained  in  this  study  agree   with  the 
previous  findings  in  other  studies  that  the  lateral  moment 
can  be  very  significant  for  ships  in  waves.   Some  definite 
trends  of  the  moment  responses  can  be  observed  from  figures 
III.l  to  III. 4.   From  the  results  in  short-crested  seas 
shown  in  figures  III. 2,  3  and  4,  one  sees  that  horizontal 
and  vertical  moments  change  very  little  with  wave  direction 
in  moderate  seas.   On  the  other  hand,  there  is  a  broad  peak 
away  from  the  head  and  following  seas  for  all  moment  responses 
in  low  sea  states.   In  severe  seas,  the  maximum  response  for 
M^  and  M   are  in  the  head  and  following  seas  while  M   stays 
almost  constant  regardless  of  wave  heading.   For  the  long 
crested  seas  of  figure  III.l,  one  can  see  that  the  vertical 
moment  peaks  at  about  60   for  low  to  moderate  sea  states  and 
assumes  extreme  values  at  head  and  following  seas  in  high 
seas.   The  equivalent  combined  moment  has  two  peaks  in  all 
cases.   The  peaks  occur  in  wave  headings  near  60   and  120 
with  a  shift  away  from  the  beam  seas  case  as  the  sea  grows 
more  severe. 
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Vertical    wave    bending    moment 
amplitude  in  regular  bow  waves 

Figure    IV. la 


Vertical  wave  bending  moment 
amplitude  in  regular  beam  and 
quartering  waves 

Figure    IV . lb 
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Horizontal  wave  bendint*  moment 
amplitude  in  regular  bow  waves 


Figure    IV. 2a 


Horizontal  wave  bending  moment 
amplitude  in  regular  beam  and 
quartering  waves 

Figure    IV. 2b 
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One  observes  from  figure  III. 2  and  III. 3  that  the 
horizontal  moments  exceed  the  vertical  moments  in  almost  all 
headings  for  low  and  moderate  seas.   However,  in  severe  seas, 
for  which  structural  design  is  based,  one  sees  the  customary 
relation  of  the  vertical  moment  exceeding  the  horizontal 
moment,  especially  in  head  and  following  seas.   The  reason 
for  this  phenomenon  is  that  lateral  bending  is  more  respon- 
sive to  higher  frequency  excitation  than  vertical  bending. 
One  can  see  from  figures  IV. 1  and  IV. 2  that  the  spectral  peak 
of  the  lateral  moment  response  for  all  wave  headings  are  gen- 
erally situated  at  higher  frequencies  than  the  corresponding 
response  curves  for  the  vertical  moments.   One  can  note  fur- 
ther that  the  spectral  ordinates  of  the  RAO  for  wave  headings 
70  ,  110   and  130   are  higher  for  the  lateral  moment  than  the 
vertical  moment.   Therefore,  when  the  energy  of  the  exciting 
waves  are  concentrated  in  the  higher  frequency  range  (i.e., 
in  low  sea  states) ,  the  lateral  moment  exceed  the  vertical 
moment.   In  high  sea  states,  the  reverse  trend  is  true. 

For  a  similar  reason  the  peaks  of  the  response  curves 
for  the  vertical  moment  occur  away  from  head  and  following 
seas  in  low  sea  states.   The  shift  in  the  location  of  the 
spectral  peak  of  the  RAO's  in  the  direction  of  higher  cir- 
cular frequency,  w,  can  be  explained  by  considering  the 
effective  wavelength  in  oblique  seas.   It  is  generally  known 
that  the  wavelength  nearly  equal  to  the  ship  length  causes 
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Figure    IV . 3 
The    Effective    Wavelength    in    Oblique    Waves 
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the  maximum  midship  bending  moments.   In  oblique  seas  the 
ship  senses  the  effective  wavelength,  L  ,  as  shown  in 
figure  IV. 3.   The  relation  with  the  actual  wavelength  is 


L   =  L/cos  £ 


where  L  =  effective  wavelength 


L  =  actual  wavelength 

f>  =   heading  angle 

As  p   approaches  90   from  either  0   or  180  ,  the  actual 
wavelength  decreases  while  the  value  of  L   stays  fixed. 
Therefore  the  peak  frequency  of  the  response  goes  up,  thereby 
shifting  the  RAO  spectrum  nearer  the  sea  spectral  energy 
density  peak  in  low  sea  states. 

The  theoretical  calculations  predict  greater  vertical 
moment  response  in  following  seas  than  in  head  seas,  contrary 
to  experimental  results.   Wahab's  results  [15],  shown  as 
figure  IV. 4,  confirm  this  apparent  discrepancy  between  the 
theory  and  experiments. 

Finally,  one  notes  that  the  combined  moment  in 
"realistic" seas  is  significantly  higher  than  either  of  the 
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components  in  all  cases  but,  because  of  the  relative  phase 
angles,  is  not  equal  to  the  sum  of  their  individual  peaks. 
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IV. 2   Comparison  of  the  Values  of  the  Combined  Moment 

Table  III.l  shows  the  comparison  of  the  response 
RMS  value  of  the  combined  moment,  M  ,  as  obtained  by  the  two 
methods  described  in  chapter  two.   It  can  be  readily  seen 
that  the  results  are  not  significantly  different,  the  aver- 
age for  all  sea  spectra  being  on  the  order  of  eight  per  cent 
(8%) .   This  confirms  the  statement  that  the  average  value  of 
the  correlation  coefficient,  /^   ,  already  considers  the  phase 
difference  of  the  vertical  and  horizontal  bending  moment 
peaks.   As  would  be  expected,  the  third  set  of  values  for  M 
were  closer  to  the  computer  values.   It  is  realized  that  the 
two  values  of  /?   that  were  used  were  derived  in  a  different 
manner  and  perhaps  should  not  be  directly  compared.   In  de- 
riving the  second  value,  it  was  found  that  the  c^^.__J-cient 
varies  strongly  with  the  state  of  the  sea  and  the  wave  di- 
rection and  therefore  that  no  single  value  will  be  "true"  for 
all  conditions  at  sea.   If  this  method  is  to  be  adopted, 
there  is  need  to  study  the  correlation  coefficient  further. 
There  is  a  dearth  of  published  material  on  the  subject. 

The  calculated  value  of  the  combined  moment  also  de- 
pends on  the  relative  phase  angle  of  the  vertical  and  hori- 
zontal moments.   Therefore,  it  is  necessary  that  the  pre- 
dicted values  be  at  least  as  accurate  as  the  predictions  for 
the  vertical  and  horizontal  moment  responses.   Comparison  of 
theoretical  and  experimental  phase  angles  shows  that  the 
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agreement  is  in  general  quite  good  [22],   The  maximum  dis- 
crepancies occurred  for  wave  headings  of  50   and  110   for 

the  wave  lengths  A/L  =  1.2  and  -  =  1.5.  (Figures  IV. 6  and  IV. 7) 

Jj 


In  view  of  the  above  mentioned  considerations  it 
would  seem  to  be  indicated  that  both  method  1  and  method  2 
to  obtain  the  combined  moment  are  equally  "good".   The  first 
method  is  probably  more  convenient  for  use  in  analytical 
studies. 
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Wave  Bending  Moment  at  Midship 
for  Series  60,  Block  0.80 
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Phase  Angle  Between  Horizontal  and  Vertical 
Wave  Bending  Moment  at  Midship 
for  Series  60,  Block  0.80 
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IV. 3   Comparison  Study  of  the  Vertical  Bending  Stress 

Data  on  the  full-scale  measurement  of  vertical 
bending  stresses  on  the  deck  of  the  UNIVERSE  IRELAND  can  be 
found  in  [25].   A  comparison  between  calculated  and  measured 
values  was  done  by  plotting  the  average  RMS  peak-to-peak 
amidship  bending  stresses.   The  analytical  values  were  ob- 
tained by  calculating  the  response  to  each  of  nine  sea 
spectra  corresponding  to  the  observed  periods  in  the  tabu- 
lated data  of  [24  3.   The  significant  wave  heights  are  the 
wind-to-wave  conversion  shown  in  table  III. 3.   The  average 
RMS  peak-to-peak  response  is  the  weighted  mean  of  the  cal- 
culated response  to  each  sea  specified  by  the  significant 
height.   This  is  best  described  by  the  relation 


?     h 

mean  RMS  =  \      Ey 2  •  p (T) 

i=l 


where 


Ev   =  RMS  peak-to-peak  bending  stress 


p(T)  =  relative  frequency  of  occurrence  of  each 
observed  period  for  a  given  significant 
wave  height. 
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The  wind-wave  relation  as  described  by  the  Pierson- 
Moskowitz  random  and  ITTC  1966  line,  the  corrected  Roll  curve 
for  the  North  Atlantic  and  the  corrected  Yamanouchi  curve  for 
the  North  Pacific  were  examined.   The  conversion  is  summar- 
ized in  table  III. 3. 

The  use  of  the  Pierson-Moskowitz  Random  line  pro- 
duced the  worst  case  in  the  comparison  study.   While  the 
wind-to-wave  relation  is  intended  to  convert  wind  measure- 
ments to  significant  wave  heights  for  seas  in  random  stages 
of  partial  development,  there  is  a  bias  in  favor  of  seas  that 
are  near  full  development  [25].   Consequently,  the  values 
based  on  this  wind-to-wave  conversion,  coupled  with  the  known 
tendency  of  the  theory  to  over-predict  [21],  are  iu^cn  higher 
than  the  measured  stresses.   The  Roll  and  Yamanouchi  curves 
were  based  on  observed  wave  heights  and  corrected  for  signif- 
icant heights.   They  are  believed  to  be  better  representation 
of  an  actual  sea.   Since  the  ship  route  is  neither  entirely 
in  the  North  Atlantic  nor  in  the  North  Pacific,  one  sees  fairly 
good  agreement  with  results  based  on  both  the  Roll  and 
Yamanouchi  relations  for  the  lower  half  of  the  sea  states  as 
described  by  the  beaufort  scale.   It  can  be  seen  that  theoret- 
ical results  over-predicted  the  average  values  of  the  RMS 
stress,  specially  in  the  higher  sea  states.   There  are  several 
probable  causes  of  this: 
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1.  Since  the  largest  divergence  of  compared  values  are 
consistently  in  the  more  severe  seas,  one  is  led 

to  believe  that  the  assumption  of  linear  motions 
for  the  ship,  in  particular  resonance  in  roll, 
exaggerates  the  analytical  values.   Obviously, 
heave  and  pitch  are  large  in  heavy  seas.   The 
divergence  caused  by  roll  resonance  gets  bigger 
as  the  energy  of  the  waves  are  concentrated  more 
and  more  in  the  lower  frequencies.   Faltinsen  [2  2] 
used  a  first  approximation  of  roll  at  resonance 
as  input  values.   It  is  also  to  be  noted  that 
the  ship  was  not  provided  with  bilge  keels  in 
the  computer  analysis.   It  is  not  known  if  the 
acutal  ship  has  these  dampeners. 

2.  Loukakis  [21]  suggested  that  deflections  of  the 
hull  can  relieve  some  bending  stresses,  so  that 
even  if  bending  moment  predictions  are  accurate 
the  resulting  conversion  into  stresses  may  result 
in  over-prediction. 

3.  The  inexact  weight  load  distribution  used  in  the 
analysis  could  have  caused  some  of  the  over  predic- 
tion.  A  qualitative  measure  of  this  effect  can  be 
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gained  by  the  excess  in  value  in  the  calmer  seas 
when  the  contribution  of  inertia  effects  to  the 
dynamic  loads  are  relatively  greater. 

4.  Exaggeration  of  the  reported  occurrence  of  the 
higher  period  waves.   In  high  seas,  when  the  ship 
is  pitching  and  heaving,  it  is  understandable 
that  such  errors  could  be  made  by  untrained  ob- 
servers.  In  this  study,  this  effect  would  be 
very  minor. 

5.  Inaccuracy  in  the  measured  data  due  to  insuffi- 
cient number  of  records  in  extreme  weather. 

For  direct  comparison  of  the  scatter  of  the  derived 
values  of  the  stress  with  the  measured  values,  it  was  necess- 
ary to  determine  the  standard  deviations  with  respect  to 
wind  measurements.   One  can  see  from  table  III. 4  that  the 
calculated  response  to  the  nine  sea  spectra  with  the  same 
significant  wave  height  were  more  widely  scattered  than  the 
measured  RMS  peak-to-peak,  values. 
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IV. 4   Bending  Moment  Trends  and  Wave  Load  Prediction 

One  can  observe  from  figures  III. 10  and  III. 11  that 
the  combined  moment  in  short-crested  seas  are  considerably 
higher  than  the  vertical  moment  in  head  seas.   The  difference 
in  severe  seas  is  in  the  order  of  15-20  per  cent.   Therefore, 
at  least  for  this  ship,  calculations  based  on  vertical  bend- 
ing only  could  lead  to  disastrous  results. 

Comparison  of  the  horizontal  and  vertical  bending 
moments  in  short-crested  seas  have  been  partially  discussed 
in  section  IV. 1.   It  was  explained  that  in  low-to-moderate 
seas  the  lateral  bending  response  to  wave  excitation  is 
greater  than  vertical  bending.   One  consequence  of  this  is 
its  possible  significance  in  fatigue  failures.   Another  is 
that  over  any  length  of  time  the  expected  wave  bending  load 
will  be  greater  for  lateral  moments  than  for  vertical  moments. 
This  is  so  because  the  ship  will  encounter  more  low-to-moderate 
seas  than  high  seas  in  its  lifetime.   The  result  for  the 
short-term  analysis  for  the  North  Atlantic  confirms  this 
statement.   The  statistics  reported  in  [24]  were  collected 
for  a  period  of  seven  years. 

Table  III. 5  shows  that  in  short-crested  seas  the 
bending  moment  varies  only  slightly  with  speed.   This  state- 
ment holds  true  for  all  the  moment  responses,  with  M^. 
having  the  greatest  variation  for  all  three  sea  states  shown. 
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This  indicates  that  reducing  speed  in  a  seaway  reduces 
dynamic  stresses  due  to  slamming,  etc.  but  will  not 
significantly  reduce  the  primary  stresses  due  to  wave 
bending  loads. 


Chapter  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  subject  of  the  combined  effect  of  horizontal  and 
vertical  wave  moments  has  been  investigated.   A  number  of 
conclusions  can  be  drawn  from  the  results  of  the  study. 

1.   The  results  obtained  by  the  MIT  Seakeeping 

Program  agreed  well  with  published  results  of 
theoretical  and  experimental  investigations 
of  ship  motions  and  loads  in  a  seaway.   The 
following  bending  moment  trends  were  -^served: 

i.   The  horizontal  moments  experienced  by  the  UNIVERSE 
IRELAND  are  not  negligible  in  comparison  with  the 
vertical  moments  in  almost   all  conditions  at  sea. 
In  general,  the  relative  significance  of  the  lat- 
eral moments  increase  with  the  size  of  the  ships 
[15].   In  low  sea  states,  the  lateral  moment 
exceed  vertical  moments . 

ii.   The  variation  with  wave  headings  of  the  moment 

response  seem  to  depend  on  the  state  of  the  sea. 

For  long-crested  seas,  the  maximum  vertical 

moment  occurs  at  a  wave  heading  near  60   for  low 

to  moderate  seas  and  at  0   and  180   in  severe 
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seas.   The  horizontal  and  combined  moments,  on 
the  other  hand,  exhibit  two  peaks  near  60   and 
120   in  all  sea  conditions.   For  short-crested 
waves,  all  moment  responses  show  a  broad  peak 
away  from  0   and  180   in  low  seas.   The  vertical 
and  horizontal  moments  flatten  out  in  moderate 
seas,  and  the  vertical  and  combined  moment 
assume  their  biggest  values  at  0   and  180   in 
severe  seas. 


The  above  behavior  of  the  moments  with  respect  to 
wave  headings  can  be  explained  by  the  character  of 
the  RAO. 

2.   The  values  of  the  combined  moment  obtained  by  the 
two  methods  described  in  this  paper  were  not  very 
different.   Therefore,  both  methods  can  be  con- 
sidered valid.   The  first  method  involves  develop- 
ing the  RAO  for  the  equivalent  combined  moment  by 
a  strip-theory  based  computational  method  and  per- 
forming the  usual  power  spectral  analysis.   Invest- 
igations on  the  correlation  between  analytical  and 
experimental  values  have  led  to  the  general  accep- 
tance of  the  method  first  introduced  by  St.  Denis 
and  Pierson.   However,  it  is  felt  that  further 
investigations  and  refinements  are  necessary  be- 
fore the  analytically  derived  values  can  be  used 
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in  actual  design.   The  second  method  makes  use 
of  the  correlation  coefficient,  /^   ,  to  take  into 
account  the  phase  difference  of  the  vertical  and 
horizontal  moments.   There  is  also  further  need 
to  investigate  this  important  statistical 
parameter. 

3.  The  comparison  study  of  the  predicted  and  measured 
vertical  bending  stresses  showed  good  agreement 

in  the  lower  half  of  the  weather  spectra  en- 
countered in  service.   In  higher  seas,  the  excess- 
ive predicted  values  are  believed  to  be  caused  by 
non-linear  motions,  contrary  to  the  assumption  of 
the  theory.   No  doubt  the  addition  oi  bilge  keels 
in  the  analysis  of  the  ship  would  have  reduced  the 
errors. 

For  the  ship's  route,  both  the  Roll  and  Yamanouchi 
corrected  wind-to-wave  relations  seemed  to  be  valid 
representations  of  the  seaway.   In  view  of  the  ex- 
pected higher-than-actual  values  by  theory,  Roll's 
curve  is  deemed  the  better  approximation  of  the 
service  route  if  used  in  conjuction  with  the  tabu- 
lated data  of  [24],  as  was  done  in  this  study. 

4.  It  was  found  that  there  was  a  significant  diff- 
erence between  the  values  of  the  vertical  moment 
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response  in  long-crested  head  seas  and  the 
combined  moment  response  in  short-crested  seas. 
Again,  it  is  stated  that  this  may  not  be  true 
for  other  ships,  particularly  the  smaller  ones. 
For  large  tankers  with  the  general  character- 
istics of  the  UNIVERSE  IRELAND,  the  contribution 
of  the  horizontal  moment  to  the  longitudinal 
stresses  must  be  considered  in  the  calculation 
of  the  primary  strength. 

Speed  is  not  a  significant  factor  in  the  magnitude 
of  the  wave  loads  acting  on  ships. 

It  remains  to  be  said  that  the  absolute  values  of  the 
combined  moment,  or  the  total  longitudinal  stress,  as  pre- 
dicted by  the  methods  in  this  study  could  not  be  compared 
with  experimental  results.   However,  the  moment  trends  were 
studied  and  found  to  be  reasonable  and,  on  the  whole,  in 
conformity  with  other  researchers1  results.   Therefore,  one 
has  confidence  on  these  derived  values.   The  theoretical  vert- 
ical bending  stress  predictions  could  have  been  closer  to 
measurements  but,  again,  one  has  firm  ideas  about  the  pro- 
bable reasons  for  the  divergence  and  therefore  one  is  not 
completely  dissatisfied. 

The  procedure  followed  in  this  study  is  fairly  in- 
volved and  the  cost  is  not  cheap.   For  further  analysis 
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of  ships  in  related  studies,  it  is  recommended  that  the  RAO, 
or  an  approximation,  be  obtained  in  the  cheapest  manner. 
Hopefully,  with  some  anticipation  of  the  expected  results, 
more  realistic  values  can  be  obtained. 
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Appendix  A 


SAMPLE  COMPUTER  OUTPUT  DATA 


I.   Irregular  Wave  Results  -  Long  Crested  Unidirectional  Seas 
Ship  Speed  =  12  kts.    Heading  Angle  =  22.5  degrees 
Significant  Wave  Height  =  19.2  feet 
Peak  Spectral  Frequency  =  0.5197 

Motions  at  Origin: 

RMS  H(l/3) H(l/10) H(l/1000) 


Heave 

0.5974 

2.3897 

3.0469 

4.6003 

Pitch 

.2456 

0.9822 

1.2524 

1.8908 

Sway 

.2137 

0.8549 

1.0901 

1.6458 

Roll 

.4386 

1.7546 

2.2371 

3.3775 

Yaw 

.1301 

0.5204 

0.6635 

1.0017 

Dynamic  Loadings : 

Vertical  B.M, 


Sta 


RMS 


H(l/3) 


H(l/10) 


H(l/1000) 


1 

0.0 

0.0 

2 

51.1 

204.4 

3 

1454.4 

5817.9 

4 

8894.8 

35579.2 

5 

41643.6 

166574.5 

6 

89928.7 

359715.0 

7 

144749.5 

578998.0 

8 

199801.3 

799205.5 

9 

251274.3 

1005097.2 

10 

296087.0 

1184348.0 

11 

329147-6 

1316590.0 

12 

347970.3 

1391881.0 

13 

353323.6 

1413294.0 

14 

344392.0 

1377568.0 

15 

321096.7 

1284387.0 

16 

284297.1 

1137188.0 

17 

180498.0 

721992.2 

18 

65363.6 

261454.4 

0 
7 
9 
5 

6 


0 

260 

7417 

45363 

212382 

458636.6 

738222.5 

1018987.0 

1281499.0 

1510043.0 

1678653.0 

1774648.0 

1801950.0 

1756399 

1637593 

1449915 

920540 

333354 


0.0 

393.6 

11199.6 

68490.0 

320656 

692451 

1114571 

1538470 

1934812 

2279869 

2534437 

2679371 

2720591 

2651818 

2472444 

2189088 

1389835 

503299 
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Vertical  B.M.  (cont.) 
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Sta. 


RMS 


H(l/3) 


H(l/10) 


H(l/1000) 


19 
20 
21 
22 


24263.0 
3266.2 
3990.7 
3683.9 


97052.3 
13065.0 
15963.1 
14735.9 


123741, 
16657 
20353 
18788, 


186825.7 
25150.2 
30729.0 
28366.7 


Horizontal  B.M. 


Sta. 

RMS 

H(l/3) 

H(l/10) 

H(l/1000) 

1 

0.0 

0.0 

0.0 

0.0 

2 

4568.3 

18273.5 

23298.8 

35176.6 

3 

15469.4 

61877.6 

78894.0 

119114.4 

4 

24915.5 

99662.0 

127069.1 

191849.4 

5 

52746.4 

210985.7 

269006.8 

406147.5 

6 

89010.3 

356041.5 

453952.9 

685379.8 

7 

129957.7 

519831.0 

662784.5 

1000674.6 

8 

170987.7 

683951.0 

872037.5 

1316605.0 

9 

208736.5 

834946.2 

1064556.0 

1.0 

10 

238483.0 

953932.0 

1216263.0 

,9.0 

11 

258256.3 

1033025.5 

1317107.0 

■  1988574.0 

12 

268522.5 

1074090.0 

1369464.0 

2067623.0 

13 

269648.6 

1078594.0 

1375208.0 

2076294.0 

14 

260647.3 

1042589.5 

1329301.0 

2006984.0 

15 

242403.2 

969613.0 

1236256.0 

1866504.0 

16 

214430.1 

857720.5 

1093593.0 

1651111.0 

17 

132947.2 

531789.0 

678031.0 

1023693.7 

18 

48485.1 

193940.3 

247274.0 

373335.2 

19 

22369.9 

89479.8 

114086.7 

172248.6 

20 

6595.0 

26380.2 

33634.7 

50781.8 

21 

456.7 

1826.8 

2329.2 

3516.7 

22 

118.8 

475.5 

606.2 

915.3 

Combined  B.M. 


Sta. 

RMS 

H(l/3) 

H(l/10) 

H(l/1000) 

1 
2 
3 
4 
5 

0.0 

4177.5 

15031.4 

30139.6 

86181.8 

0.0 

16710.1 

60125.7 

120558.4 

344727.5 

0.0 

21305.4 

76660. 3 

153712.0 

439527.5 

0.0 

32167.0 

115742.1 

232075.0 

663600.3 

Combined   B.M.     (cont.) 
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Sta. 

RMS 

H(l/3) 

H(l/10) 

H(l/1000) 

6 

164735.4 

658941.7 

840150.7 

1268462.0 

7 

253106.8 

1012427.2 

1290844.0 

1948922.0 

8 

342644.4 

1370577.0 

1747486.0 

2638362.0 

9 

424373.1 

1697492.0 

2164303.0 

3267672.0 

10 

492267.6 

1969070.0 

2510565.0 

3790460.0 

11 

540433.3 

2161733.0 

2756210.0 

4161336.0 

12 

566890.0 

2267560.0 

2891139.0 

4365052.0 

13 

571413.7 

2285655.0 

2914210.0 

4399885.0 

14 

554399.9 

2213199.0 

2821829.0 

4260409.0 

15 

513669.5 

2054678.0 

2619714.0 

3955255.0 

16 

452582.7 

1810331.0 

2308172.0 

3484887.0 

17 

281527.8 

1126111.0 

1435792.0 

2167764.0 

18 

99230.1 

396920.7 

506073.9 

764072.3 

19 

37856.3 

151425.4 

193067.4 

291494.0 

20 

6591.3 

27805.4 

35451.8 

53525.4 

21 

4085.7 

16342.9 

20837.2 

31460.1 

22 

3671.7 

14687.0 

18725.9 

28272.5 
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II.   Irregular  Wave  Results  -  Short-crested  Multi-Directional 
Seas 

Ship  Speed  =  12  kts . 

Significant  Wave  Height  =  4  3.3  feet 

Peak  Spectral  Frequency  =  0.3461 

Principal  Wind  Direction  =  22.5  degrees 

Motions  at  Origin: 

RMS         H(l/3)       H(l/10)   H(l/1000) 


Heave 

6.2544 

25.0175 

31.8974 

48.1588 

Pitch 

1.1184 

4.4734 

5.7036 

8.6113 

Sway 

3.5923 

14.3690 

18.3205 

27.6603 

Roll 

4.1180 

16.4720 

21.0018 

31.7086 

Yaw 

0.6277 

2.5109 

3.2014 

4.8334 

Dynamic 

Loadings : 
Vertical  B. 

,M. 

Sta 

RMS 

H(l/3) 

H(l/10) 

•  H(l/1000) 

9 

660410.1 

2641640.0 

3368092. 

.0 

5085158.0 

10 

770714.8 

3082859.0 

3930646. 

.0 

5934504.0 

11 

852958.4 

3411833.0 

4350088, 

.0 

6567779.0 

12 

902672.4 

3610689.0 

4603629. 

.0 

6950577.0 

13 

916354.1 

3665416.0 

4673406. 

.0 

7055927.0 

14 

891059.9 

3564239.0 

4544406. 

.0 

6861161.0 

15 

827132.9 

3308531.0 

4219378. 

.0 

6368923.0 

16 

728960.0 

2915840.0 

3717696. 

,0 

5612991.0 

17 

456989.5 

1827958.0 

2330646. 

.0 

3518819.0 

Horizontal  B.M. 


Sta 


RMS 


H(l/3) 


H(l/10) 


H(l/1000) 


9 
10 
11 
12 
13 
14 
15 
16 
17 


440864. 

518902, 

575644. 

606.50, 

609593, 

584602, 

535193, 

463863 

273484 


1763456.0 

2075610.0 

2302577.0 

2424600 

2438373 

2338410 

2140775 

1855452 

1093937 


2248406, 
2646403, 
2935786, 
3091365, 
3108925, 
2981473, 
2729488, 
2365701, 
1394770, 


3394653, 
3995549, 
4432461, 
4667354, 
4693867, 
4501440, 
4120991, 
3571745, 
2105830, 
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Combined    B.M. 


Sta 

RMS 

H(l/3) 

Hd/10) 

H(l/1000) 

9 

958354. 

,0 

3833416.0 

4887606. 

.0 

7379326.0 

10 

1121474. 

.0 

4485896.0 

5719517. 

,0 

8635349.0 

11 

1242414. 

.0 

4969656.0 

6336311. 

.0 

9566587.0 

12 

1310208. 

.0 

5240832.0 

6682061. 

,0 

10088601.0 

13 

1320890. 

.0 

5283560.0 

6736539. 

.0 

10170852.0 

14 

1273295. 

.0 

5093180.0 

6493804, 

,0 

9804371.0 

15 

1170845. 

.0 

4683380.0 

5971309. 

.0 

9015506.0 

16 

1020343. 

.5 

4081374.0 

5203752. 

.0 

7856645.0 

17 

616199, 

,8 

2464799.0 

3142619, 

,0 

4744738.0 

Appendix  B 


CALCULATIONS  FOR  AN  AVERAGE  VALUE  OF  pr 


VH 


For  comparison  with  the  value  of  the  correlation 
coefficient  of  /0   =  0.32  as  obtained  from  figure  III. 5, 
an  average  of  /0VH  was  calculated  from  the  moment  responses 
computed  by  equation  (6)  of  the  first  method. 


M 


T 


-  E 


V 


-  K2E 


M 


H 


VH 


2K  E 


V 


M 


H 


k 


As  the  long-term  extreme  wave  moments  were  not  ex- 
trapolated, the  mean  Attt  value  was  obtained  by  averaging 
over  wind/wave  directions  and  for  representative  sea  states 
The  results  of  the  calculation  are  summarized  below: 
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V  =  16  knots 

Short-crested  seas  responses 


90 


wind/wave  direction 

h=7.4  ft. 
/VH 

h=24.3  ft. 
1   VH 

h=43.3  ft. 
1    VH 

o°  n 

22.5° 
45.° 
67.5° 
90°  o 

112.5 

1350 

157.5° 

180° 

0.771 
.799 
.633 
.743 
.676 
.587 
.460 
.293 
.214 

.692 
.700 
.697 
.650 
.555 
.441 
.384 
.414 
.468 

.590 
.598 
.603 

.367 
.334 
.368 
.395 

average  ^VR 

0.575 

0.556 

0.477 

mean    for    low,    moderate   and   high    seas 
TVH   =    0.536 


Appendix  C 

CALCULATION  OF  THE  MEAN  RMS  PEAK-TO-PEAK 
VERTICAL  BENDING  STRESS  AND  THE 
STANDARD  DEVIATION  WITH  RESPECT  TO  WIND 

The  calculated  bending  moment  was  converted  into  stress 
by  the  relation: 

Stress  =  Moment  *  Section  Modulus 

where  the  section  modulus  of  the  ship  at  the  deck  is  given  as 
566,794  in2  -  ft  [253. 

The  computer  output  gave  the  RMS  of  the  response.   To 
get  the  peak-to-peak  RMS  value,  the  data  was  multiplied  by 
the  factor  2/2. 

The  mean  value  of  the  response  to  a  sea  described  by 
its  significant  height  is  given  by: 


9 

my  =  I    V.  p(T.) 


i=l 


where  tcl       =  mean  RMS  of  the  peak-to-peak  response  to  a  sea 
specified  by  the  significant  height. 


V.  =  RMS  of  the  peak-to-peak  response  to  the  ith 
sea  spectra 
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p^T.^  =  probability  of  occurrence  of  the  ith  sea  spectra. 

p(_T.)  is  approximated  by  the  relative  frequency  of 
occurrence  of  the  different  observed  wave  periods  for  a 
given  wave  height  in  the  tabulated  data  of  [24]. 

Average  values  for  all  the  areas  in  the  ship's  route 
were  obtained.   The  calculation  is  illustrated  in  table  A.l 
for  the  one  case  of  h  =  13.2  feet  (wave  height  code  05  in 
[24])  . 

The  standard  deviation  of  the  vertical  bending  stress 
with  respect  to  wave  measurements  is  calculated  from  the 
relation 


S22  "  J^i  *  %f   P(Ti) 


where  S2  =  standard  deviation  of  the  bending  stress  with 
respect  to  waves. 

For  direct  comparison  with  experiments   based  on 
wind  measurements,  the  following  relation  is  used: 


2      2 
Sl   =  S2   + 


fs  2  -  Ms 

[S3     12 


tan2e  [23] 
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The  numerical  calculation  is  illustrated  below  for 
the  Roll  wind- to-wave  relation. 
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Beaufort  No 

h(ft) 

m 

V 

2/      r  ^ 

s   (ton-ft) 

s   (kpsi) 

s  (kpsi) 

1 

7.2 

123079. 

3.465xl09 

0.573 

.239 

2 

7.6 

158320. 

5.616xl09 

.878 

.296 

3 

8.0 

186316. 

6.36xl09 

.0992 

.315 

4 

9.0 

229732. 

5.188xl09 

.0810 

.285 

5 

10.0 

305014. 

7.572xl09 

.1184 

.344 

6 

12.4 

443627. 

1.322xl010 

.2064 

.454 

7 

15.0 

541361. 

2.506xl010 

.391: 

5 

8 

18.2 

803314. 

3. 279x10 10 

.5120 

.716 

9 

22.0 

964491. 

5. 530x10 10 

.8637 

.929 

10 

26.0 

1156994. 

7. 499x10 l0 

1.1712 

1.082 



2     2  J  r  2   ah 

Sl  =  S2  +  [S3  "  71" 


1    2 

J  tan  G 


Beaufort  No 

h(ft) 

2 
s2  (kpsi) 

H(ft) 

Un6l  ft  j 

s  (kpsi) 

s1(kpsi) 

1 

7.2 

.0573 

.8 

.115692 

2.5 

.374 

2 

7.6 

.0878 

.4 

.115692 

2.5 

.414 

3 

8.0 

.0992 

.6 

.115692 

2.5 

.435 

4 

9.0 

.810 

1.0 

.115692 

2.8 

.430 

5 

10.0 

.1184 

2.5 

.115692 

3.4 

.516 

6 

12.4 

.2064 

2.5 

.115692 

4.3 

.668 

7 

15.0 

.3912 

3.5 

.115692 

5.1 

.758 

8 

18.2 

.5120 

3.8 

.115692 

6.1 

.997 

9 

22.0 

.8637 

4.0 

.115692 

7.0 

1.225 

10 

26.0 

1.1712 

4.0 

.115692 

7.8 

1.403 
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